Non-equilibrium vertex correction: disorder scattering in magnetic tunnel junctions 
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We report a first principles formalism and its numerical implementation for treating quantum 
transport properties of nanoelectronic devices with atomistic disorder. We developed a nonequilib- 
rium vertex correction (NVC) theory to handle the configurational average of random disorder at the 
density matrix level so that disorder effects to nonlinear and nonequilibrium quantum transport can 
be calculated from atomic first principles in a self-consistent and efficient manner. We implemented 
the NVC into a Keldysh non-equilibrium Green's function (NEGF) based density functional theory 
(DFT) and applied the NEGF-DFT-NVC formalism to magnetic tunnel junctions with interface 
roughness disorder. Our results show that disorder has dramatic effects to nonlinear spin injection 
and tunnel magneto-resistance ratio. 
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Quantitative understanding of impurity effects is cru- 
cial for nanoelectronics where device properties are 
strongly influenced by or even built on such effects. Ex- 
amples are electron scattering by dopants in semiconduc- 
tor nanowires[l| and field effect transistors, spin scat- 
tering by disorder in magnetic tunnel junctions^, and 
transport of spin polarised current in dilute magnetic 
semiconductors 0. Unintentional impurities sit inside 
a device at unpredictable locations and therefore, any 
physical quantity predicted by theory should be aver- 
aged over impurity configurations. In ab initio calcula- 
tions, one may carry out this average by generating many 
impurity configurations for a given concentration a;, cal- 
culating the relevant physical quantity for each configu- 
ration, and finally averaging the results. Such a brute 
force calculation is often not practical for at least two 
reasons. First, when x is small as is typical in the case of 
semiconductor devices, say 0.1%, one would need a thou- 
sand host atoms to accommodate just one impurity atom. 
Second, it is known that impurity average may require 
huge number of configurations These requirements 
make a calculation prohibitively large. Considerable ef- 
fort has therefore been devoted in literature to develop 
approximate techniques which avoid brute force. In this 
regard, a widely used technique is the coherent potential 
approximation (CPA)0] as implemented in KKrQ and 
LMTO0 first principles methods. So far CPA has been 
applied to equilibrium electronic structure and trans- 
port calculations 0. However, most nanoelectronic de- 
vices operate under nonequilibrium conditions, for in- 
stance one wishes to predict nonlinear current-voltage (I- 
V) characteristics. It is thus very important to develop 
appropriate nonequilibrium techniques for impurity av- 
eraging. 

Here we report our solution of the atomistic nonequilib- 
rium impurity average problem for quantum transport. 
We start from a state-of-the-art real space atomistic 



quantum transport formalism where density functional 
theory (DFT) is carried out within the Keldysh nonequi- 
librium Green's function (NEGF) framework[i[i3|. The 
basic idea of NEGF-DFT is that the device Hamiltonian 
and electronic structure are determined by DFT, the non- 
equilibrium quantum statistics of the device physics is 
determined by NEGF, and the transport boundary con- 
ditions under external bias are handled by real space nu- 
merical technique. We deal with impurity average at sin- 
gle particle retarded Green's function level by CPA[1], 
and at NEGF level by evaluating a nonequilibrium vertex 
correction (NVC) term. The NEGF-DFT-NVC formal- 
ism allows us to construct nonequilibrium density matrix 
self-consistently that includes impurity averaging. We 
then apply our NEGF-DFT-NVC formahsm to investi- 
gate effects of interface roughness disorder in a magnetic 
tunnel junction (MTJ). Our results indicate that disorder 
effect can drastically and qualitatively infiuence nonlin- 
ear I-V curves and tunnel magneto-resistance ratio. 

We consider a two-probe device consisting of a scatter- 
ing region and two semi-infinite leads extending along the 
transport direction z to z = ±oo, as shown in Fig[TJ The 
system is periodically extended along the transverse (x,y) 
direction. Note the scattering region includes several lay- 
ers of lead atoms 0. A bias voltage 14 is applied across 
the leads to drive a current flow, i. e. jii— = eVb where 
are electrochemical potentials of the left/right leads. 
We assume that impurities exist inside the scattering re- 
gion randomly but not in the leads. We further assume 
that any atomic position R in the scattering region may 
be occupied by two atomic species, the host and impurity 
atoms labelled hy Q — A,B with concentrations and 
C| such that -h C| = 1. 

In a NEGF-DFT self-consistent analysis fo^ of ordered 
systems, the non-equilibrium density matrix is calcu- 
lated by NEGF G<(£;), i.e. fi{E) - G<{E), here 
G< satisfies the Keldysh equation G< = G'^Y.<G^ 
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where = (G'^)^ is the retarded Green's function. 

= iTifi + iTrfr is the lesser self-energy where fi^r 
are Fermi functions of the left/right leads; F; are line- 
width functions describing coupling of the scattering re- 
gion to the leads and they can be calculated by standard 
iterative methods 0]. Translational invariance of ordered 
system allows one to evaluate all quantities in the unit 
cell by integration over two dimensional (2D) Brillouin 
zone (BZ) in the (x,y) direction. When there are impu- 
rities, translational symmetry is broken. The spirit of 
CPA is to construct an effective medium theory by an 
impurity configurational average that restores the trans- 
lational invariance. For NEGF this means calculating 
G = G^S<G-^. Even though there is no impurity in 
the leads to affect S^, the impurity average (• • •) cor- 
relates Green's functions with G-^. In particular, 
G ^ G S<G due to multiple scattering by the impu- 
rities. To calculate G , we introduce a self-energy FAryc 
that is a consequence of impurity scattering at nonequi- 

librium, such that G — G (E^ + ^nvc)G . T^vc 
is called nonequilibrium vertex correction (NVC) whose 
equilibrium counterpart is well known in calculations of 
Kubo formula by Feynman diagrammatic techniques 11 1 . 
There is however a major qualitative difference here: 
r NVC depends on the non-equilibrium quantum statisti- 
cal information of the device scattering region while the 
equilibrium VC does not. 

We found that T^vc is most conveniently calculable 
using a site oriented calculation scheme, for this reason 
we develop our NEGF-NVC theory within TB-LMTO 
DFT implementation [13], using CPA 0,|l]to describe the 
averaged system. In this approach, impurity average for 
any single site physical quantity Xr is given by 



Q=A,B 



(1) 



where is the conditional average over a particular 
atomic specie Q at site R which is calculated by 



Xr, = nlXRic\ 



(2) 



Here r\^ is the occupation of site R by atomic specie Q 
and its average is rf^ = G^ . Eq. ([T]) means the average of 
a quantity at site i? is a linear combination of contribu- 
tions from each atomic species. 

The technical derivation details are given in the Sup- 
plemental Material associated with this paper (13|, here 
we briefly outline the spirit of the theory. The impurity 
average of the site diagonal NEGF is carried out by appli- 
cation of Eq.ll]), i.e. G^j; — ^ ^r^rr ■ The atom 

Q=A,B 

resolved NEGF G^^ gives atom resolved average local 
charge density ~ G^f, which is needed in the DFT 
self-consistent iterationsfol]. To find G^^ , we use Eq.®. 



The final expressions of G^^ are given in Eqs. (10,28,29) 
of the Supplemental Material [l3|, and it is related to 
Tmvc of the auxihary NEGF. The calculation of NVC 
is carried out by application of single-site approximation 
(SSA) within CPA-based multiple scattering theory as 
summarised in Supplemental Material[13'!| where the fi- 
nal expression is given by Eq.(23) there. From T^vc we 
obtain G^^ hence the averaged density matrix for 
atom Q on site R. The charge density is used to cal- 
culate device Hamiltonian for the next step in the DFT 
iteration and this procedure is repeated until numerical 
convergence. 

An extremely stringent test of our NEGF-DFT-NVC 
formalism and its numerical implementation is carried 
out by calculating two-probe devices at equilibrium, and 
check if the fluctuation-dissipation relationship is satis- 
fied or not. Mathematically, fluctuation-dissipation theo- 
rem dictates G^^ — G^^ — G^'^ at equilibrium. Here, 
calculation of G^^ requires NVC while calculations of 

Gf^fj and Gj^]^ do not. For many disordered device 
structures including that in Fig[T] we have checked, the 
fluctuation-dissipation relationship is always satisfied to 
at least one part in a million and the final tiny difference 
can be attributed to numerical calculation issues. Impor- 
tantly, we found that NVC is extremely important: with- 
out it the density matrix and transmission coefficients 
can have large errors and even become qualitatively in- 
correct. 

After the NEGF-DFT-NVC self-consistent calculation 
is converged we calculate current-voltage (TV) char- 
acteristics by Landauer formula, where an additional 
vertex correction must be done on the transmission 



coefficients [16[. At low temperature the I-V curve is given 
by: 



Tr 



dE 



(3) 



The integrand of the above expression is the transmission 
coefficient where the impurity average, once again, corre- 
lates the retarded and advanced Green's functions 
{g^'^ are auxiliary Green's functions of G^ -^, see Sup- 
plemental Material) connected by line-widths Ti^r cor- 
responding to the auxiliary g. We write the averaged 
transmission into a coherent part and a vertex part which 
describes the inter-channel scattering events: 

T ^ Tr [Tig-^Trg-^] + Tr [Tig-^T'ycg^] (4) 

where T'y^ is obtained from the expression of F nvc by 
replacing with F,.. The equilibrium conductance for 
a spin channel is given by T(e^/ft-). 

As an important application of the NEGF-DFT-NVC 
formalism, we have investigated a disordered MTJ shown 
in Fig[T] which consists a vacuum (Vac) tunnel barrier 
sandwiched by two Fe leads. The Fe/Vac interface has 
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roughness disorder which substantially influences spin de- 
pendent transport at equilibrium ^1 7] . Here we focus on 
nonequilibrium. In our Fe/Vac/Fe MTJ, the left /right 
Fe/Vac interface layers has x% and (1 — x)% Fe atoms re- 
spectively, and the rest are vacuum sites. The scattering 
region consists of ten perfect atomic layers of Fe oriented 
along (100) on the left and right ending with the rough 
interface sandwiching four vacuum layers. The scatter- 
ing region is connected to perfect Fe left /right leads ex- 
tending to z = ±oo. Because the entire structure, after 
impurity average, is periodic along the transverse x, y 
directions, we found that very careful two-dimensional 
Brillouin zone (BZ) sampling is necessary in calculating 
the density matrix: we use 200 x 200 k-mesh to ensure 
excellent numerical convergence. For the I-V curve cal- 
culation, Eq.®, 300 X 300 BZ k-mesh is used for each 
point in the energy integration. For other DFT details 
we follow standard TB-LMTO literature [i3|. 

Fig[2^ is a semi-log plot of equilibrium {Vb = 0) con- 
ductance versus disorder x for both spin-up and -down 
channels, G'|,Gj^. The four curves correspond to mag- 
netic moments of the two Fe leads having parallel or 
anti-parallel configurations (PC or APC). Since the left 
interface is chosen to be FcxVaci-x while the right 
Fei^xVO'Cx^ G-\,i are both symmetric about x — 0.5 in 
PC (black squares and red circles). For APC they are 
not symmetric but satisfy G'\{x) — G|(l — x), as ex- 
pected. Impurity scattering dramatically decreases 
in PC, and has relatively less effect for G| and for APC. 
It was well known [Till that for perfect interfaces, the 
surface electronic states of Fe give resonance transmis- 
sion. These resonances are destroyed rapidly by the in- 
terface disorder as x changes from zero to 50% leading 
to the drastic reduction of G| in PC. An important de- 
vice merit for MTJ is the tunnel magneto-resistance ratio 
(TMR) defined by total tunnelling currents for PC and 
APC: TMR=(/^C' _ jAPCyjAPC_ equihbrium when 

all currents vanish, we use equilibrium conductances to 
calculate TMR. Figl^D plots equilibrium TMR versus 
X showing a dramatic effect of disorder. In particular, 
TMR drops to very small values, even to slightly neg- 
ative values, as x is increased from zero. These equi- 
librium features are consistent with previous super-cell 
calculations 17 1 . 

We now investigate nonequilibrium properties when 
Vb ^ so that current flows through. To show the impor- 
tance of NVC, we have calculated I-V curves at a; = 0.05 
by including vertex correction only at the level of trans- 
mission coefficient, i.e. without NVC in the NEGF-DFT 
self-consistent iterations of the density matrix: the solid 
lines (green) in Figl3^ plot this result. In comparison, 
the dashed lines (red) plot the full results where NVC is 
included. The substantial differences indicate that NVC 
is extremely important for obtaining correct results at 
nonequilibrium. Fig[3}D and its inset plot TMR versus 
bias for four values of x = 0.0, 0.05, 0.3, 0.5, obtained by 



the full NVC formalism. For zero or small values of x, 
TMR reduces with Vb as is often seen in experimental 
measurements [l9j| . For larger x, for instance x 0.5, 
TMR can go negative as Vb is increased. Indeed, ex- 
perimental measurements had seen[2^ negative TMR at 
large Vb, although for different MTJs and possibly differ- 
ent physical origin. Very dramatically, at a; = 0.3, the 
entire TMR curve is negative: here the absolute value 
of TMR actually increases with Vb (see inset). These 
behaviours of TMR strongly suggest that interface dis- 
order play very important roles for nonequilibrium spin 
injection. 

Fig HI plots spin currents and TMR versus disorder 
X a,t Vb = 0.544V. This is to be compared with FigH 
where Vb — 0. A finite bias breaks left-right symmetry 
of the atomic structure and therefore, the spin currents 
do not have a symmetric behaviour about a; = 0.5 any- 
more. Both spin currents (Figd^) and TMR (FigUja) 
varies with disorder x in substantial ways. In particular, 
TMR rapidly dips to negative values when x is increased 
to about 20%. So far we have focused on devices where 
the left has a Fea;Vaci_j; interface while the right has 
Fei_:rVac^. We have also apphed the NEGF-DFT-NVC 
formalism to devices where the left and right interfaces 
are disordered totally differently. The inset of Fig[4|3 
plots TMR for such a system where left interface has 
X — 0.3 while the right interface has x = 0.05. For this 
system TMR is negative and its absolute value decreases 
as Vb is increased which is qualitatively similar to what 
discussed above. 

In summary, we have developed a nonequilibrium ver- 
tex correction theory and its associated software for 
analysing quantum transport properties of disordered na- 
noelectronic devices at nonequilibrium. The impurity av- 
eraging of the nonequilibrium density matrix is facilitated 
by the NVC self-energy that is related to quantum sta- 
tistical information of the device scattering region. Our 
NEGF-DFT-NVC theory has several desired features, 
including atomistic first principle, non-equilibrium, ef- 
ficient configurational average and self-consistency. This 
allows us to analyse nonequilibrium quantum transport 
of realistic device structures including realistic atomic 
substitutional impurities. Using this tool, we have cal- 
culated nonlinear spin currents in Fe/Vac/Fe MTJ with 
interface roughness disorder, and found that effects of 
NVC can play a dominant role in determining the prop- 
erties of spin injection. 
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FIG. 1: (Colour online) Schematic of atomic structure of the 
Fe/Vac/Fe magnetic tunnel junction. The two Fe/Vac inter- 
faces have roughness disorder. Fe: yellow spheres; vacuum: 
white spheres. 




FIG. 2: (Colour online) (a) Conductance G|,| versus disorder 
X at equilibrium. Red circles: Gj in PC; black squares: G| 
in PC. Blue down-triangles: G; in APC; Green up-triangles: 
Gt in APC. (b) TMR Versus x. 
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FIG. 3: (Colour online) (a) Comparison of I-V curves with 
disorder x — 0.05. Solid lines (green): current for PC (up- 
triangles) and APC (down-triangles) without using NVC in 
density matrix self-consistent iteration. Dashed lines (red): 
current for PC (circles) and APC (squares) using the full NVC 
formalism, (b) TMR versus bias voltage Vb for four different 
values of x. The main figure is for x=0.0 (black squares) 
and x=0.05 (red circles); the inset for x=0.3 (black down- 
triangles) and x=0.5 (red up-triangles) 
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FIG. 4: (Colour online) (a) Spin currents versus disorder x 
at bias Vb = 0.544V, for PC and APC. Red circles and black 
squares: spin currents for spin-up and -down in PC; green up- 
triangles and blue down-triangles: spin currents for spin-up 
and -down in APC. (b) TMR versus x at the same 14. Inset of 
(b): TMR versus Vh for a device where left and right interfaces 
have different values of a;, on the left interface x — 0.3, on the 
right X = 0.05. 



